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SUMMARY 

An algorithm Is presented which can be used to develop compliance matrices 
for cracked bodies. The method relies on the numerical solution of singular 
Integral equations with Cauchy-type kernels and provides an efficient and accu- 
rate procedure for relating applied loadings to crack opening displacements. 

The algorithm should be of Interest to those performing repetitive calculations 
In the analysis of experimental results obtained from fracture specimens. 


INTRODUCTION 

Two dimensional problems In linear elastic fracture mechanics are often 
reduced to a singular Integral equation (or system of Integral equations) of 
the form 

f f t ± M 1 * /] K(x ' t} f(t) dt = g(x) - 1 < x < 1 (D 

where f(t) Is a function to be determined, and K(x,t) and g(x) are known 
functions related to the geometry of the cracked body and the loading on the 
crack surface(s), respectively. The function f ( t ) , which Is called the dis- 
location density. Is the slope of the crack profile, and. If the crack Is 
closed at Its tips, satisfies the consistency condition 

/] f(t) dt = 0 (2) 

In this paper an algorithm Is presented which can be used to develop a compli- 
ance matrix for a cracked solid when such a formulation Is used. This matrix 
relates the stresses at the collocation points arising from the applied loads 
(Including those applied to the crack surfaces) to the values of the opening 
of the crack at the Integration points. The method relies on the numerical 
procedure developed by Gerasoulls (ref. 1), which Is used to reduce 
equation (1) to a system of algebraic equations for unknown values of f(t) 
at discrete points In the Interval [-1,1]. 


FORMULATION 

The algorithm Is best explained through a simple example. For the cracked 
plate shown In figure 1 the governing equations for the dislocation density are 
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where y Is the shear modulus, k = 3 - 4t> for plane strain, and v Is 
Poisson's ratio. 


The crack opening displacement Is given by 


u + - u" = J f(t) dt 


(5) 


The exact solution to this problem Is 
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Following [1], (after nondlmenslonallzlng by setting 2y/ira 00 (ic + 1) equal 
to unity), f(t) Is expressed by 

f(t) = 4»(t) (1 - t 2 )" 1 / 2 (8) 

and <j>(t) Is approximated as piecewise quadratic In [-1,1]. The result Is 
that equations (3) and (4) are reduced to a system of algebraic equations 
through quadrature formulas. The details of the quadrature can be found In 
reference 1 and are omitted here. The results are 
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We note that the unknown values $ are those at the Integration points t^, 
while the stresses on the right hand side of equation (9) are at the colloca- 
tion points X|<. 

For Illustration we take five points for the quadrature, and equations (9) 
and (10) become 
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In equation (11), matrix G Is characteristic of the geometry of the prob- 
lem and matrix L represents the loading and crack closure condition. What Is 
usually of Interest in such a problem Is the stress Intensity factor, and since 
the stress Intensity factor Is proportional to the value of $ at the end- 
points (ref. 1), equation (11) Is solved for the unknown vector $. 


Instead of doing this, In this paper the Inverse of matrix G Is obtained, 
and the product of this matrix and ( 1 — 1^> -1 /2 integrated term by term 

to obtain a matrix C which Is called the compliance matrix for this particu- 
lar geometry. The Integration of each term Is performed using the weights for 
the Lagrange Interpolation polynomials, since the function <t> Is approxi- 
mated In this manner. The results of the Integration lead to the following 
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We note that matrix C Is not a square matrix. This Is because the number of 
Integration points Is one more than the number of collocation points. Matrix 
relates the displacement of the crack faces at the point t^ to the 
stress at point x^. 

For the present problem the applied loading vector Is unity, and premulti- 
plying It by the compliance matrix leads to 

<u + - u"), = 0.0 

(u + - «f) 2 = 0.276 

(u* - u~) 3 = 0.318 

(u + - u ) 4 = 0.276 

(u + - u~) 5 » 0.0 

The above displacements agree with the exact solution. 
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The usefulness of the compliance matrix becomes evident when one Is Inter- 
ested In Investigating the effects of many loading cases, and In particular, 

If the stresses along the crack surfaces depend on the crack opening displace- 
ments. Such loadings are used In models for fiber reinforced concrete, rocks, 
ceramics, and other materials where microcracking, fiber bridging, and other 
nonlinear effects are modeled as nonlinear springs along the crack surfaces 
(refs. 2 to 5). For these models equation (1) becomes nonlinear, and an Iter- 
ative solution Is needed. With the use of the compliance matrix, this Iterat- 
ion procedure Is efficient and fast. 


As an example, let us assume that the crack surfaces are bridged by fibers, 
and that the stresses transmitted by the fibers to the crack depend on the 
opening of the crack. The displacements along the crack will be governed by 
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(14) 


where now the stresses are decomposed Into those arising from the applied 
loads, and those due to the fiber bridging. The function <j(u + - u~) Is 
determined from experiments (ref. 6). For the first Iteration the stresses due 
to fiber bridging are assumed to be equal to zero. Premultiplying the stresses 
arising from the applied loading by the compliance matrix results In the first 
approximation to the crack opening displacements. From these displacements the 
first approximation to the fiber bridging stresses are determined, and these 
are applied to the crack surfaces. The procedure Is repeated until convergence 
Is reached. This procedure was used In reference 4 (where experiments per- 
formed on concrete and fiber reinforced concrete were analyzed) and convergence 
was observed to be very fast (only several Iterations were needed for three and 
four point bending specimens). 


CONCLUSIONS 

An algorithm has been presented which can be used to develop compliance 
matrices for cracked bodies. The usefulness of the matrices becomes evident 
when one Is Interested In performing parameter studies to Investigate the 
effects of various loadings on crack opening displacements. Even though the 
example presented In this paper Involved Mode-I loading, the method can be 
extended to Include mixed mode problems (Including three-dimensional problems). 


REFERENCES 

1. Gerasoulls, A.: The Use of Piecewise Quadratic Polynomials for the 

Solution of Singular Integral Equations of Cauchy Type. Comput. Math. 
Appl . , vol. 8, no. 1, 1982, pp. 15-22. 

2. Keer, L.M.; and Meade, K.P.: A Note on a No-Slip Interface Crack, 

J. Appl. Mech., vol. 49, no. 2, June 1982, pp. 454-455. 

3. Wecharatana, M.; and Shah, S.P.: A Model for Predicting Fracture 

Resistance of Fiber Reinforced Concrete. Cem. Concr. Res., vol. 13, 
no. 6, Nov. 1983, pp. 819-829. 


4 



4. Ballarlnl, R.; Shah, S.P.; and Keer, L.H.: Crack Growth In Cement Based 

Composites, Engineering Fracture Mechanics, vol . 20, no. 3, 1984, 

pp. 433-445. 

5. Catalano, D.M.; and Ingraffea, A.R.: Concrete Fracture: A Linear Elastic 

Fracture Mechanics Approach, Report No. 82-1, Dept, of Structural 
Engineering, Cornell University, 1982. 

6. Shah, S.P.; Stroeven, P.; Dalhulsen, D.; van Stekelenberg, P.: Complete 

Stress-Strain Curves for Steel Fiber Reinforced Concrete In Uniaxial 
Tension and Compression. Proc. Int. Symp., RILEM-ACI-ASTM, Sheffield, 
1978, pp. 339-408. 


5 





1. Report No. 

NASA CR-1 79478 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Compliance Matrices for Cracked Bodies 

S. Report Date 

July 1986 


6. Performing Organization Code 

7. Authors) 

R. Ballarlnl 


6. Performing Organization Report No. 

None 


10. Work Unit No. 

9. Performing Organization Name and Address 



Cleveland State University 
Department of Civil Engineering 

11. Contract or Grant No. 

NCC 3-46 

neveiana, umo id 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

National Aeronautics and Spj 
Washington, D.C. 20546 


Contractor Report 

ice Administration 

14. Sponsoring Agency Code 

506-42-11 ( E— 31 58) 

IS. Supplementary Notes 

Final Report. Project Manager, Anthony M. Calomlno, Structures Division, NASA 
Lewis Research Center, Cleveland, Ohio 44135. 


16. Abstract 


An algorithm Is presented which can be used to develop compliance matrices for 
cracked bodies. The method relies on the numerical solution of singular Integral 
equations with Cauchy-type kernels and provides an efficient and accurate proce- 
dure for relating applied loadings to crack opening displacements. The algorithm 
should be of Interest to those performing repetitive calculations In the analysis 
of experimental results obtained from fracture specimens. 


17. Key Words (Suggested by Authors)) 

Integral equations 
Compliance matrix 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 39 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of pages 

22. Price* 


*For sale by the National Technical Information Service, Springfield, Virginia 22161 























National Aeronautics and second class mail 

Space Administration 

Lewis Research Center address correction requested 

Cleveland, Ohio 44135 

Official Business 
Penalty (or Private Use $300 



Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





